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Phase diagrams are the most reliable way of predicting phase changes in a system. In this work, we create the 
phase diagram of the Ga( 20 _ x ) Al x + bimetallic cluster system, for which the monometallic clusters display greater-than- 
bulk melting behaviour. Employing first-principles Born-Oppenheimer molecular dynamics in the microcanonical 
ensemble, we describe the solid-liquid-like phase transition in two distinct compositions GanAlg + and GasAli 7 + . The 
clusters show peaks in specific heat at 824 K and 922 K respectively which is above the corresponding bulk alloy 
melting temperatures. Mean squared displacements, diffusion coefficients and atoms-in-molecules analysis emphasise 
the differences between the environments of the internal and surface atoms. An Arrhenius-like description of the 
‘ease’ with which internal and surface sites (or atoms) can exchange positions with temperature is used to quantify 
the greater-than-bulk melting character which corroborates the view of a Ga( 20 _ x )Al x + cluster being a remnant of the 
corresponding bulk alloy material with additional characteristics specific to its size and composition. 


1 Introduction 

Nanoalloys comprising two or more metal atoms within 
a finite framework offer novel physical and chemical 
properties significantly different from their bulk ahoy 
counterparts^]. The sensitivity of their properties to sys¬ 
tem size and composition gives ample room for fine tuning 
as desired for various applications in catalysis, optics and 
electronics[22]. Understanding the effects of size, struc¬ 
ture, composition, chemical ordering and temperature on 
the emergent behaviour of nanoalloys thus plays a key role 
in their fabrication and functional optimization. 

Although a lot of scientific literature within current 
nanoscience addresses the role of ground state nanoalloy 
configurations on the observed properties p~8l l22t f2| H7] . ef¬ 
fects at finite-temperature have received less attention [T3] 
HU- One such property is the solid-liquid-like phase tran¬ 
sition temperature which has already been found, both 
experimentally [20, [21] and in theoretical simulations [31, 
06] [471 [45], to depend on the size and composition of 
nanoalloys. Moreover, the phase transition behaviour 
changes qualitatively as well; different melting stages (pre¬ 
melting) were observed in Cu-Au system[32] in contrast to 
Ag-Au bimetallic clusters [33] which showed a single stage 
melting. Complex structural transformations, including 
switching among homotops before complete melting, have 
also been observed!9, 1 . 

Compositional dependence of the phase transition tem¬ 
perature in alloys is studied through the construc¬ 
tion of phase diagrams. Extension of the phase di¬ 


agram to nanoscale materials for predicting the solid- 
liquid-like phase transition temperature has been already 
attempted pSi 51, 39 23]. Models typically used for bulk 
phases [T2l [30] are unsuitable for these systems and hence, 
sensitive first-principles treatment of the electronic struc¬ 
ture is necessary. In our previous study [57 ], we showed 
that density functional theory (DFT) was able to cap¬ 
ture the solid-liquid-like phase transition temperature of 
GaigAl + in excellent agreement with experiments [35. 
In this work, we explore two more cluster sizes in the 
Ga( 20 _ x )Al x + cluster series: GanAlg + and Ga 3 Ali 7 + , and 
attempt to create a phase diagram of Ga-Al nanoalloy 
system using the phase transition temperature data of 
monometallic phases, Ga 2 o + and Al 2 o + respectively [50l 

EH]. 

Cationic gallium clusters have been observed to show 
melting-like transition at temperatures significantly above 
the bulk melting temperature of gallium [TUI 29, [11] 9i T9]. 
However, a recent study by Jarrold and co-workers[43] 
have found the gallium cluster comprising 95 atoms to 
be the first cluster in the size range 60-183 atoms to melt 
below the bulk melting temperature. Although not ob¬ 
served in experiments due to limitations on the temper¬ 
ature range 36 , greater-than-bulk melting behaviour has 
also been predicted for aluminium clusters [38]. An alu¬ 
minium atom, being isoelectronic and having a similar 
ionization potential (IP), has a vacant d-shell and is 5.4% 
smaller in atomic radius in comparison to a Ga atom. 
Moreover, the overlap of valence s and p orbitals is differ¬ 
ent for both the metals since the 4 s band of Ga is shifted 
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to a higher binding energy in comparison to the 3s band of 
AIDS]. These differences raise the question as to how the 
thermodynamics of Ga-Al clusters change with composi¬ 
tion in the nano-regime, and also how is it different from 
the corresponding bulk Ga-Al alloy system. The aim of 
this report is to answer the above question and also to un¬ 
derstand further the differences between the behaviour of 
the internal and surface atoms of the cluster and how they 
affect the phase transition, as has been observed in Ga 2 o + 
and Al 2 o + clusters [38]. The organisation of the work is 
as follows: section 2 describes the computational method 
adopted for this study. In section 3 we present the results 
obtained using various statistical tools and indices to ex¬ 
plain the spatial and temporal variations occurring during 
the phase transition of GanAlg” and GasAl^y clusters and 
go on to complete the phase diagram of the 20-atom Ga- 
Al system. Finally we summarise and conclude in section 
4. 


2 Computational details 

We adopted a similar methodology for studying the melt¬ 
ing transitions in GanAlg + and GasAli 7 + clusters as has 
been described in our previous works|37l 138] . Density 
functional theory (DFT) calculations were performed us¬ 
ing the plane waves basis set based code Vienna ab-initio 
simulation package (VASP)[26j 27, 54] |25]. Large core 
(i.e. only 3s 2 Sp 1 and 4 s<2 4p 1 electrons were considered as 
valence for A1 and Ga respectively) projector-augmented 
wave (PAW) [8j [28] pseudo-potentials in combination with 
the generalized gradient approximation in the Perdew- 
Wang form (GGA-PW91) [51] 42, were used. The start¬ 
ing structures for GanAlg + and GasAli 7 + clusters (as 
shown in Fig. [5]) were determined using the stacked plane 
(SP)[50l[37] Ga 2 o structure. Each of the 20 available sites 
in the SP configuration of Ga 2 o was substituted by a sin¬ 
gle A1 atom and based on the calculated free energy of the 
cluster, the sites were numbered in ascending order. The 
first nine and seventeen sites served as the positions for A1 
atoms in GanAlg + and Ga 3 Aln + respectively. 

Born-Oppenheimer molecular dynamics (BOMD) as im¬ 
plemented in VASP in addition to parallel tempering (PT) 
was used to perform the molecular dynamics simulations 
in the microcanonical ensemble. 19 different temperatures 
between 250 K and 1150 K at 50 K intervals were cho¬ 
sen to capture the solid-liquid-like phase transition. Rela¬ 
tive convergences of the canonical and microcanonical spe¬ 
cific heat curves obtained using the multiple histogram 
(MH) method [T4l 5(1, ascertained using the sliding win¬ 
dow analysis [50 , were achieved after a total run time of 
160 ps (step size 2 fs) for GanAlg + and 46.2 ps (step 


size 0.8 fs) for GasAli 7 + . A smaller step size was used 
in GasAli 7 + to counter the increase in energy during mi¬ 
crocanonical runs. Parallel tempering (PT) helps to avoid 
the dependence of thermodynamics on the starting struc¬ 
ture and the possibility of the structure getting trapped 
in one of the local minimas in the potential energy sur¬ 
faces of GanAlg + and Ga 3 Aln + . After every 100 time 
steps, each of the 19 different microcanonical runs was 
stopped. Two configurations were randomly selected and 
based on the calculated acceptance probability [15], swap¬ 
ping was either accepted or rejected. The spatial charge 
distribution analysis was performed using freely available 
software [3] based on the quantum theory of atoms-in- 
molecules (QTAIM)jJ. 


3 Results and Discussion 

3.1 Phase transition in GanAlg and 
Ga 3 Al+ 

The top panel of Fig. |T] shows the canonical heat capacity 
curves obtained from the MH method for GanAlg + and 
GasAli 7 + clusters normalized to the classical specific heat 
capacity (Co = (3 N — 6 + 3/2 )ks where N is the number 
of atoms and ks is the Boltzmann constant) accounting 
for rotational and vibrational contributions. The solid- 
liquid-like phase transition temperature, corresponding to 
the peak in canonical specific heat capacity plot, is 824 K 
for GanAlg + and 922 K for Ga 3 Aln + cluster respectively. 

Analyses of the indices pertaining to the cluster struc¬ 
ture at different temperatures help to relate the struc¬ 
tural changes occurring in the cluster to its thermody¬ 
namic behaviour. Melting, in the context of clusters, is 
a change from the solid-like state to the liquid-like state. 
The root mean squared (RMS) bond-length fluctuation 
(Berry parameter[52]) identifies this phase transition i.e. 
the melting-like temperature as the temperature above 
which the average change in bond-lengths converge. Ex¬ 
pressed as S rrns , 




(< r ij >t~< Tij >?) 1/2 


where N is the total number of atoms, is distance be¬ 
tween atoms i and j and <> t denotes the time average, 
this index of fluctuation quantifies the response of bond 
vibrations to increasing energy (or temperature). 

Shown in the bottom panel of Fig. [l]are the calculated 
firms values for GanAlg + (shown in blue) and Ga3Aln + 
(shown in red) clusters. The fluctuations in bond-length 
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Figure 1: Comparison of canonical specific heat curves of 
GanAlg + and GasAli 7 + (top panel) and the root-mean- 
squared (RMS) bond-length fluctuations (bottom panel) 
as a function of temperature. 

are quite low (~ 0.1 — 0.15) at the lowest energy and start 
to increase with temperature. However, above the cal¬ 
culated solid-liquid-like phase transition temperature of 
824 K in GanAlg + and 922 K in GasAli 7 + , the values 
converge to (^ 0.3) indicating a structural transition. It 
must, however, be stated that although at low and high 
temperatures the values of 5 rms are consistent with PT 
independent simulations, the absolute value of the index 
in the intermediate temperature range is subject to the 
influence of PT. 

Mean squared displacement (MSD) helps to compare 
the relative mobilities of the internal and surface atoms. 
Calculated as per Eq. [2j 

1 M 

< r i{t) > = — ^i(^ra)] 2 (2) 

m= 1 


Lowest Energy 



t(ps) 


Figure 2: Comparison of mean squared displacement 
(MSD) of GanAlg + at lowest (246 K) and highest en¬ 
ergy (931 K) as a function of time. Atoms in red and blue 
are the central A1 atoms; green and orange corresponds 
to MSD averaged over all the A1 atoms in the central and 
top ring respectively; indigo represents the average MSD 
of Ga atoms as also shown in inset. For references in the 
text, the more highly coordinated internal atom (red) is 
referred as a and the other (coloured blue) as b. 


where M is the number of time origins, the individual 
atomic contributions to the overall behavior of the sys¬ 
tem can be analyzed. Fig. [2] shows the MSD of gallium 
and aluminium atoms of GanAlg + at the lowest and high¬ 
est energy. Coloured in indigo, the amplitude of the MSD 
of the surface gallium atoms (averaged) does not change 
significantly as the energy is increased. A1 atoms present 
at the internal sites (red and blue) show the maximum in¬ 
crease in the amplitude of MSD from lowest to the highest 
energy. Moreover, at the highest energy all the A1 atoms 
have almost similar MSDs indicating a liquid-like phase. 
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Figure 3: Comparison of mean squared displacement 
(MSD) of GasAli 7 + at lowest (286 K) and highest energy 
(1079 K) as a function of time. Atoms in red and blue are 
the central A1 atoms; turquoise and indigo corresponds to 
MSD averaged over all the surface A1 and bottom ring Ga 
atoms respectively. Orange represents the MSD of top Ga 
atom. For references in the text, the more highly coordi¬ 
nated internal atom (red) is referred as a and the other 
(coloured blue) as b. 


Fig. [3] shows the MSD calculated for GasAli 7 + at the 
lowest and highest energy. In contrast to the Ga atoms in 
Gan Alg + (average MSD of all Ga atoms coloured indigo in 
Fig. [2]) where the MSD have similar amplitudes for all the 
atoms at the lowest and highest energies, one of the three 
Ga atoms situated at the top (colored orange in inset) has 
MSD values between those of internal atoms and other 
surface atoms. This difference, however, vanishes as the 
energy is increased and all the surface atoms have similar 
MSD values. It is for the internal A1 atoms (red and blue) 
that the MSD changes with increasing temperature. At 



Figure 4: Diffusion coefficient of internal and surface 
atoms as a function of temperature. The vertical dotted 
line denotes the melting temperature. 


the highest energy, all the atoms (Ga and Al) have MSD 
values reflecting a liquid-like state in GasAli 7 + . 

The diffusion coefficient helps to quantify the mobility of 
atoms in the cluster. It is calculated as per Eq. [3] where D 
is the diffusion coefficient of the specific atom, < r* (t) > is 
the mean squared displacement of the corresponding atom 
calculated as per Eq. [2] and d is the dimensionality of the 
cluster system (which is 3). For small systems such as 
clusters, the diffusion coefficient is meaningful only over 
a small time range where the MSD achieves a constant 
slope [7]- 


D = _L lim ±i±M± 

2 d t—>-oo t 


(3) 


Shown in Fig.0]are diffusion coefficient values calculated 
for GanAlg + and GasAli 7 + clusters in the top and bot¬ 
tom panels respectively. The MSD values over which the 
slope is calculated have been kept fixed between 1 and 2 
A 2 . MSD values of all the surface gallium and aluminium 
atoms have been averaged as shown in black and green 
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for both the clusters. In both cases, the diffusion coeffi¬ 
cients of the surface atoms are very high with the trend 
increasing with temperature. The average coordination 
number, calculated using the first minimum in pair distri¬ 
bution function g(r) as cutoff distance for coordination, of 
GanAlg + (Fig. 1, Supplementary Information) shows that 
the internal atoms (red and blue) occupy an internal site 
at all temperatures below the melting temperature of 824 
K. In the top panel of Fig. 0J the diffusion coefficient of 
both the internal atoms is much lower at all temperatures 
up to 701 K after which the internal atom a (coloured 
red in inset of bottom panel, Fig. El shows fluctuations 
in GanAlg + . However, it is only after the melting tem¬ 
perature that the diffusion coefficient of the internal atoms 
shoots up and at the highest average temperature becomes 
comparable with the rest of the A1 atoms in the Gan Alg + 
cluster. 

Although the surface Ga and A1 atoms of Ga 3 Aln + 
show similar behaviour with increasing temperature, con¬ 
stant exchange among atoms to occupy an internal site at 
temperatures closer to the melting temperature (as also 
clear from the average coordination number plot in Fig. 2, 
Supplementary Information) makes it much harder to tag 
a particular atom as internal. To counter this situation, we 
attribute the term ‘ persistent internaV to an atom which 
satisfies the following two conditions: (a) the average co¬ 
ordination number has to be greater-than-or-equal-to 9.5 
and, (b) the time during which condition (a) is satisfied 
during the continuous MD run has to be more than 1.5 
ps. The second condition also helps to make sure that 
an atom occupying an internal site at a particular energy 
run is favourable rather than being an artefact of parallel 
tempering. 

Coloured orange in the top and bottom panels of Fig.0] 
the diffusion coefficient of the ‘ persistent internaV atom(s) 
in both the clusters is much lower in comparison to the 
surface atoms below the melting temperature. Above the 
melting temperature, none of the atoms in GanAlg + sat¬ 
isfy the condition to be considered as c persistent internaV. 
In GasAli 7 + , in contrast, the ‘ persistent internaV atom 
criteria is fulfilled up to 967 K, above the melting tem¬ 
perature of 922 K. Moreover, in GanAlg + the diffusion 
coefficient of the ‘ persistent internaV atom is an average 
of both the internal atoms (red and blue) up to 582 K. 

The Quantum Theory of Atoms in Molecules 
(QTAIM)[4| aims to define atomic boundaries within a 
molecule using the topology of the total electronic charge 
density of the molecular system. The volume occupied 
by an atom, which is an open system, is governed by the 
topological condition of zero-flux of the electronic charge 
density (vp( r ) ' n i r ) = 0)- Integrating the electron den¬ 
sity within each atomic volume gives the total electronic 



(a) Ga+ (b) Gai 9 Al+ (c) GanAl 9 + 


(d) Ga 3 Ali 7 + (e) Al 20 + 

Figure 5: Starting structures in the capped sphere (GaJo) 
and stacked plane (GaigAl + , GanAlg + , Ga 3 Aln + and 
Al 2 o + ) configurations with A1 atoms in white and Ga 
atoms in cyan. 

charge associated with the corresponding ion. 

Shown in Fig. 3 (Supplementary Information) are the 
partial charges associated with each of the gallium and 
aluminium atoms in GanAlg + and Ga 3 Ali 7 + respectively. 
There is a strong charge segregation between the internal 
and surface atoms in both the cases. Moreover, among 
the charged surface atoms it is the gallium atom which 
becomes negatively charged and aluminium which be¬ 
comes positively charged in both cases. A similar picture, 
where an electrostatic cage confines the internal atoms 
was also observed for the monometallic Ga 2 o + and Al 2 o + 
clusters |38'. 

3.2 Phase diagram of Ga-Al 

Incorporating the results from our previous 
investigations[37, [38] of Ga 2 o + , Al 2 o + and GaigAl + , 
we have created a binary phase diagram of the 20 atom 
Ga-Al nanocluster system. Irrespective of the bimetallic 
cluster’s composition, the solid-liquid-like phase transition 
temperature is always above the corresponding bulk alloy 
melting temperature [M] as also shown in Table [TJ 

The energy of mixing, calculated as per Eq. 0] re¬ 
flects the energy gain (or loss) for a mixed cluster relative 
to the pure clusters. A negative value indicates mixing 
i.e. formation of the nanoalloy is energetically favourable, 
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Figure 6: Excess energy for all the compositions in 
Ga 2 o-ccAl cc + , panel (a), and the corresponding atoms-in- 
molecules analysis where the red dot, blue open circle and 
turquoise circle represents negative, neutral and positive 
charges. 


whereas a positive value indicates segregation between the 
participating atoms. 


Table 1: Comparison of the melting temperatures of Ga- 


A1 alloys in cluster and bulk phase (M]. 



Cluster T m (K) 

Bulk T m (K) 

AT m (K) 


616 

303 

313 

Gai 9 Al+ 

686 

353 

333 

GanAlg + 

824 

646 

178 

GasAli7 + 

922 

857 

65 

ai 20 + 

993 

933 

60 


771 _ j-,Ga(2o — x )Al+ X -pAl^Q 20 X Ga^ 0 

X^exc — Xj tot 20 20 r " t ° t 

Fig. [6] (a) shows the energy of mixing for all the cluster 
compositions in the series Ga^o-z) Al cc + . Adding an alu¬ 
minium atom to Ga^ 0 gives a positive value to the energy 
of mixing indicating segregation of the A1 atom in the gal¬ 
lium dominated Gaig Al + cluster. The energy of mixing re¬ 
mains positive for GaigAl 2 + and from Gai 7 Al 3 + onwards 
it becomes negative, achieving a minimum at the most 
favourable composition, i.e. GanAlg + . This trend is also 
visible from the change in melting temperature between 
the Ga( 2 o-^)Ala, + clusters and the corresponding bulk al¬ 
loy counterparts as shown in Tabl^H The change in melt¬ 
ing temperature first increases from Ga 2 o + to Gaig Al + but 
thereafter decreases with increasing A1 composition thus 
reflecting how the energy of mixing affects the greater- 
than-bulk melting behaviour of Ga( 2 o-x)Al x + clusters. 

In order to understand further how the charge distribu¬ 
tion among atoms, i.e. with respect to internal & exter¬ 
nal sites, affects the segregation or mixing behaviour in 
each cluster composition in Ga^o-z) Al^A, we performed 
atoms-in-molecules analysis as shown in panel (b) and (c) 
of Fig. [6l The optimized structure for each of the nanoal¬ 
loy clusters was obtained from the SP configuration of 
Ga 2 o as has been explained in Section 2. From the ob¬ 
tained partial charges on each atom (<^), an atom is con¬ 
sidered neutral if —O.le < qi < O.le, positively charged if 
qi > O.le and negatively charged if qi < —O.le. The GaJo 
capped sphere structure biff has one internal negatively 
charged atom and one positively charged surface atom 
with the remaining atoms being neutral. In GaigAl + , the 
internal A1 and Ga atom becomes strongly positively and 
negatively charged respectively, thus forming a dipole with 
no negatively charged surface atoms. In GaigAl 2 + and 
Gai 7 Al 3 + , both the internal atoms are positively charged. 
However, for Gai 7 Al 3 + , there is a negative charge on the 
surface. Starting from Gai6AU + , there is one neutral A1 
atom occupying the internal site with the other being posi¬ 
tively charged until Gai 4 Al 6 + and negatively charged from 
Gai 3 Al 7 + to Gai 7 Al 3 + . Both the internal atoms are neg¬ 
atively charged for Ga 2 Alig + , GaAlig + and Al 2 o + - Corre¬ 
lating with the energy of mixing, as shown in Fig.[6^a), we 
see that segregation (a positive energy of mixing) corre¬ 
sponds to either the formation of a strong dipole between 
the internal atoms (as in GaigAl + ) or the absence of neg¬ 
ative charge at the centre or surface (as in GaigAl 2 + ). 
Moreover, for Gai 3 Al 7 + , Gai 2 Alg + and GaAlig + , along 
with the monometallic Ga 2 o + and Al 2 o + clusters, there is 
no negatively charged surface atom. Strong charge seg¬ 
regation between the internal and surface sites plays an 
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important role in the greater-than-bulk melting behaviour 
as also reflected by the diffusion coefficients (Fig. [4]) of the 
representative clusters [38]. 

The average coordination number, diffusion coefficient 
and atoms-in-molecules analysis bring out the differences 
between the internal and surface sites in these clusters. 
The internal atoms being confined in an electrostatic cage, 
breaking of which leads to the solid-liquid-like phase tran¬ 
sition in Ga( 2 o-x)Al x + clusters, is a physical picture which 
can help explain the greater-than-bulk melting charac¬ 
ter (GB) of this system. Greater-than-bulk melting be¬ 
haviour, observed in gallium 9 and tin [35] clusters, is an 
exception to the melting point depression with size as pre¬ 
dicted by Pawlow[30]. Models typically applied to un¬ 
derstand the size dependence of the phase transition in 
clusters [23, 3D have expressed the cluster melting tem¬ 
perature as a function of the bulk melting temperature 
with a parameter catering for the size dependence compo¬ 
nent (where the nanoparticle is generally assumed spher¬ 
ical). Although this may be a reasonable approximation 
for cluster sizes having hundreds of atoms, a plethora of 
research currently addresses the different spatial geome¬ 
tries assumed by the putative global minimum structures 
of nanoparticles having fewer than tens of atoms[5] ITS] . 

How different is a gallium-aluminium nanoalloy cluster 
from the corresponding bulk alloy material? Taking into 
account the GB character, energy of mixing and the as¬ 
sociated intrinsic elemental character (to be common in 
bulk and cluster), we have attempted to provide a sim¬ 
plistic view of the Ga( 20 _ x )Al x + cluster system. Our con¬ 
tention is that a cluster can be viewed as a piece of the 
bulk matter with additional properties specific to a par¬ 
ticular cluster size, as expressed in Eq. [5] for our specific 
case. Our intention is two-fold: firstly, to find a way to 
‘quantify’ the inherent and size dependent GB character 
among Ga( 20 _ x ) Al x + clusters and, secondly to compare the 
obtained melting temperature from Eq. [5] with those ob¬ 
tained from our DFT calculations to gauge the applicabil¬ 
ity of this picture. 

Tm(Ga( 2 o- x )Al+) = Tb u ik(Ga : Al :: (20 — x) : x) 

+ Vgb + Vexc (5) 

Here T m and Tbuik are the cluster and bulk alloy melting 
temperatures for similar composition and vgb, v exc are 
the functions representing the contributions of the greater- 
than-bulk melting character and the energy of mixing re¬ 
spectively. 

To ‘quantify’ the GB character, we revert back to the 
picture of the internal atoms which are confined in an 
electrostatic cage by the surface atoms in Ga( 20 _ x )Al x + 
clusters, disruption of which is required for melting. We 


view this difference in the environments of the internal and 
surface atoms as an ‘energy barrier’ which exists between 
the internal and surface sites. With increasing tempera¬ 
ture, the internal and surface atoms have enough energy to 
cross this existing ‘energy barrier’ and thus either replace 
each other or independently become a surface or internal 
atom which subsequently leads to a change in their re¬ 
spective mobilities in their new sites below the melting 
temperature. This change in the cluster either due to a 
change in the number of internal atoms or the type of in¬ 
ternal atom (considering swapping of atoms due to parallel 
tempering or an exchange of positions between the inter¬ 
nal and surface atoms) is, from here onwards, termed as 
‘ Change-Swap (CS)’. Characterising the liquid-like state 
as one where all the atoms have an equal likelihood of 
occupying an internal or surface site, we have attempted 
to understand the strength of the ‘energy barrier’ among 
the five representative structures in Ga^o-^) AG+ series 
through an Arrhenius-like equation, Eq. [6j which is also 
used to describe diffusion. 

In kcs — ^ A ^barrier (tJ^) (6) 

Here, kcs is the rate constant in ns -1 , A is the pre¬ 
exponential factor, T is the temperature in Kelvin and 
Ebarrier characterizes the energy barrier between the in¬ 
ternal and surface sites. The ratio of total number of times 
a CS occurs to the total MD trajectory length would give 
the rate of exchange (kcs) at each temperature. This al¬ 
lows us a way of using Eb arr i er as a quantifiable ‘measure’ 
of the strength between the core (internal atoms) and sur¬ 
face atoms i.e. an intrinsic, cluster-dependent and inten¬ 
sive quantity which could be correlated with the greater- 
than-bulk melting character (vgb)- 

Shown in the inset of Fig. [8] are the data points and the 
corresponding linear least squares fit based on the ‘per¬ 
sistent internal atom criterion (a)’. We observe stronger 
fluctuations in the natural logarithm of the CS rate below 
the melting temperature and thus decided to perform the 
least squares fit to data points above the melting temper¬ 
ature. The obtained Eb arr i er values from the fit has been 
tabulated in Table [2] Although we do observe a clear in¬ 
crease in the energy barrier with increasing Al atoms in the 
Ga( 20 _^)Al ;E + cluster series, we found the ‘persistent in¬ 
ternal’ atom criterion (a) quite loose for the temperatures 
below the melting temperature for tagging an atom to be 
internal The high mobility of the surface atoms led to 
frequent expansion and shrinkage of the cluster volume as 
observed in the MD. Thus, although when observed over 
a time range as in criterion (b), the ‘persistent internal’ 
atom does extremely well to tag a particular atom as in¬ 
ternal, it needs further refining when monitoring changes 
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Figure 7: Rate of change-swap (CS) between the internal 
and external atoms in the representative Ga^o-cc) AC+ 
clusters. The temperatures at which the rate of CS was 
calculated is written above the bars. 

at every time step. Furthermore, from MD observations 
we found two ways in which the internal and surface atoms 
and sites were exchanged: (i) stage-wise where a surface 
atom became internal in subsequent MD steps, and (ii) 
due to parallel tempering. With increasing temperatures, 
we also observed changes in the number of internal and 
surface atoms reflecting presence of isomeric structures 
near the melting temperatures. We detail our observa¬ 
tions from the MD movies in section 3 (Supplementary 
Information). 

To counter the problem of accurately tagging a partic¬ 
ular atom as internal at each MD step, we employed the 
convex hull technique [d which for a finite set of points (in 
this case 20) is the union of all convex combinations of 
them. For each configuration of 20 points, we perform the 
convex hull in both two and three dimensions. In three di¬ 
mensions, we disregard the atoms which form the vertices 
of the hull. In the two dimensions, we calculate the con¬ 
vex hull for the top view, front view and side view of the 
cluster. A cartesian coordinate (atomic site) is considered 
internal if it is a non-hull vertex in all the three respective 
two dimensional convex hulls along with being a non-hull 
vertex in the corresponding three dimensional convex hull 
also satisfying the criterion (a) of the ‘persistent internal^ 
atom. 

Fig 0 shows the rate of CS for the representative clusters 
during the microcanonical simulations, at average temper¬ 
atures below the melting temperature, above the melting 
temperature and at the highest energy simulated in which 


Figure 8: Arrhenius plot showing the dependence of nat¬ 
ural logarithm of the rate of CS between the internal and 
surface sites to the inverse of temperature. The dashed 
lines are the least-square straight line fits to the data 
points obtained from the convex hull method. The in¬ 
set shows the corresponding data points obtained from 
the persistent internal atom criterion (a). The line fit has 
been done only for temperatures above the melting point 
for each cluster. 

case all clusters are fully liquid-like. We have ignored those 
CS brought about by parallel tempering where there is a 
change only in the number of internal sites (isomeric) as 
opposed to the atom occupying the internal site as well 
(homotopic). In all cases, we see an increase in the rate of 
CS with temperature. However, we also observe that be¬ 
low the melting temperature the rate of CS is dependent 
on the amount of mixing among representative clusters. 

Shown in Fig. [8] is the corresponding natural logarithm 
of the rate of CS calculated at average temperatures for 
different microcanonical trajectories. The temperatures 
for which CS was less than 10 discrete instances during 
the MD trajectory, or the contribution of PT to CS was 


Table 2: Comparison of the Eb arr i er values obtained using 
the two methods described above. _ 


‘Persistent internal’ 
atom criterion (a) 

Convex hull 
method 

GaJ 0 

729.5 

2165.6 

Gai 9 Al+ 

1306.5 

3301.4 

GanAl 9 + 

1433.9 

3767.3 

Ga 3 Ali 7 + 

2570.2 

3921.3 

> 

to 

o 

+ 

2889.5 

3982.3 
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Figure 9: Comparison of the DFT calculated cluster melt¬ 
ing temperature (red) of the five representative composi¬ 
tions of Ga( 2 o-x) Al x + cluster system to that obtained con¬ 
sidering the picture described in Eq. 0 


more than 15% have been discarded, leaving a minimum 
of eight data points. Shown as dashed lines are the corre¬ 
sponding linear least squares line fit to the values. Com¬ 
paring to the inset of Fig. 0 we observe that the spread 
in data at temperatures below the melting temperature 
has been significantly reduced, reflecting the importance 
of added screening to gauge a CS situation. The slopes 
of the line-fit, i.e. Eb arr i er for the different clusters, have 
been tabulated in Table 0 We find that with increasing A1 
content in the cluster, the energy barrier for a CS situation 
increases. It must however be noted that the actual mag¬ 
nitude of Eb arr i er depends not only on the number of data 
points and thus the fitting itself but also the method used 
to calculate the rate of CS. From Table 0 correlating with 
the AT m values tabulated in Table [lj we observe that there 
is an inverse relationship between the relative temperature 
difference between the cluster and the corresponding bulk 
Ga-Al alloy and the greater-than-bulk melting character 
(GB). 

Expressing vqb = &f(Eb arr i er ), we first perform a lin¬ 
ear least squares fit to the E\y arr i er values obtained using 
the convex hull method for the pure phases, and follow by 
optimising a. Represented as green triangles we show in 
Fig. [9] the obtained melting temperatures of the five rep¬ 
resentative clusters when GB character has been added to 
the bulk melting temperatures. Furthermore, we express 
u exc = (3E exc and thereby correct the melting tempera¬ 
tures of the mixed phases which are now within 2% error 
from the DFT calculated (red circles) melting tempera¬ 
tures, shown as blue triangles in Fig. [9j 


First principles calculations predict greater-than-bulk 
melting behaviour for the mixed gallium-aluminium 
bimetallic clusters: GanAlg + and GasAli 7 + . This be¬ 
haviour has already been observed experimentally for 
Ga 2 o + & GaigAl + and predicted for Al 2 o + - Moreover, 
the melting temperature of the Ga( 2 o- x )Al x + nanoalloys 
increases with increasing A1 content. The differences be¬ 
tween the environments of internal and surface sites of 
the mixed clusters has been convincingly demonstrated, a 
feature which seems common to the Ga( 2 o- x )Al x + cluster 
series. An Arrhenius-like analysis describes this difference 
as an energy barrier which gives a way to quantify the 
greater-than-bulk melting character of these nanoalloys. 
Furthermore, we propose a picture where the Ga( 2 o-x) Al x + 
clusters are viewed as components of the corresponding 
bulk alloy phase with a systematic increase in the melt¬ 
ing temperature which can be ascribed to this finite size 
specific energy barrier. 
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